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The methods of geometrical optics are extended so that they may be applied to gyrotropic bodies. 
Various internal and external reflections are considered at nonparallel planar interfaces and means 
of determining the ray path or direction of energy flow are derived. Nonplanar geometries may be 
represented by the tangent planes at the various points of incidence. A method is given for computing 
the phases of the various fields. These values may be used to determine the reflected fields from 
such a gyrotropic body. 



1. Introduction 

A number of investigators have applied geometrical optics to anisotropic media, notably in 
ionospheric research [Budden, 1961] where ray tracing procedures have been carried out for aniso- 
tropic interfaces, and also of more recent studies where sources have been included f Arbel, 1960; 
Arbel and Felsen, 1963; Felsen, 1963a, and b] and scatterers embedded in anisotropic regions 
[Felsen, 1963a, b, 1964a, b; Rulf and Felsen, 1964]. While these studies have often involved spe- 
cial orientations of the gyrotropic axis (parallel or perpendicular to the interface) or special medium 
parameters (uniaxial case) to permit the derivation of explicit results for various radiation and dif- 
fraction problems, the relevance of these techniques to more general situations has been indicated. 

A modified geometrical optics method also has been developed for approximating the electro- 
magnetic scattering properties of isotropic bodies. This method has been successfully applied 
to a number of different cases including cyhnders, spheres, and a prolate spheroid [Peters and 
Thomas, 1962; Thomas, 1962; Kawano and Peters, 1963; Kouyoumjian, Peters, and Thomas, 1963; 
Kawano and Peters, 1964]. These papers on the modified geometrical optics method have shown 
the manner in which diffraction of rays and caustics may be treated. In order to apply these same 
methods to gyrotropic media, geometrical optics for gyrotropic bodies is required. The purpose of 
this paper is to outline a geometrical optics technique for gyrotropic bodies. Much of the material 
is well known but has not been presented in concise form such that the modified geometrical optics 
method can be readily appHed. 

Wait [1961] has developed a boundary value solution for the infinite circular gyrotropic cylinder 
with an axial rtiagnetic field. This is a special case in which the index of refraction is independent 
of the direction of propagation through the medium. The modified geometrical optics method 
[Lee, Peters, and Walter, 1964a, and b] has also been apphed to this particular case and the results 
obtained are in excellent agreement with those obtained from the boundary value solution. 

In all cases the static magnetic field is chosen so that it lies in the plane of incidence, i.e., the 
x-z plane (see fig. 1). Then this plane contains all refracted ray paths. If the static magnetic 
field had a y component, this would no longer be true. This does not represent a hmitation on the 
methods but it does reduce the complications involved in obtaining numerical results. 



' The work reported in this paper was supported in part by Contract Number AF 19(604)-7270 between Air Force Cambridge Research Laboratories, Office of 
Aerospace Research, Bedford, Mass. and The Ohio State University Research Foundation. 
2 Presently employed by Bell Telephone Laboratories, Inc., Murray Hill, N.J. 
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Figure 1. Plane wave incident on a plasma gyrotropic 
interface. 
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2. Ray Optics for Arbitrary Gyrotropic Boundaries 

The fundamentals needed to trace a ray path through a finite, homogeneous, gyrotropic body 
and to compute the magnitude and phase of the emergent rays are developed in this paper by 
considering only planar interfaces. This might appear to be a severe restriction; however, modi- 
fied geometrical optics as applied to isotropic bodies is also based on the planar interfaces. Snell's 
law and the various reflection and transmission coefficients are all derived under the assumption 
of infinite planar interfaces. Yet this modified geometrical optics method has been remarkably 
successful in computing the radar cross section for isotropic spheres and cylinders. This modified 
geometrical optics method also yields remarkably accurate radar cross sections for the plasma 
cylinder with an axial static magnetic field [Wait, 1961; Lee, Peters, and Walter, 1964 a, and b]. 
The reflection and transmission coefficients for the planar gyrotropic interfaces have been devel- 
oped previously and are given in the appendix. 

The same assumption made in developing the modified geometrical optics method in these 
previous cases is also made here for this general gyrotropic case. That assumption is that the 
planar interfaces represent the tangent planes of any curved body at the point of incidence and the 
point of reflection. The ray technique for finite homogeneous gyrotropic bodies will now be 
summarized. 

A ray is incident at the origin as shown in figure 1. The value of S is obtained from the angle 
of incidence, ^/, as S = sin Oi and C = cos di, where the incident plane wave is given by 



{// = f/^,e-jWx+c^)^ 



(1) 



The parameters of the gyrotropic plasma medium are given by X and Y which are defined by 



x=^= 






(2) 



and 



y _ CO// _ fipeHo 



0) 



mo) 



Y{Ut-\-lsk\ 



(3) 



where 



coh is the cyclotron frequency, 
o)N is the plasma frequency, 
6o is the permittivity in free space. 



(4) 
(5) 
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and 



fjio is the permeability in free space, 
A^ is the number of electrons per unit volume, 
e is the charge of electron, 
m is the mass of electron, 

Hq is the applied static magnetic field, 

/i, /'i, /.} are the direction cosines of the static magnetic field. 



i, k are the unit vectors for the coordinate system of figure 1. 



These represent the physical parameters of the gyrotropic body to be treated. If that body has 
a curved surface it should be noted that S is a function of the point at which the ray enters the body 
for a particular plane wave incidence, or conversely that the coordinate system of figure 1 is always 
chosen so that the z-axis is normal to the surface at the point of incidence which is the origin. 
Thus the values of S, C, /i, h, and /a are all functions of this angle of incidence, di. 

Once Si is found, values of S, C, /i, /2, and 4 are readily obtained. These parameters uniquely 
determine the values of the coefficients a, /3, y, 8, and e of Booker's quartic equation, which is 
given by 

F{q) = aq'' + ^q^ + yq^ ^bq-V e = ^, (6) 

where the coefficients are functions of direction, St. For the lossless case, the collision frequency 
is negligible, and 



a = {\-Y')^X{Yl-\) = [l-Y')+X{l'^''-\\ (7) 

P = 2XSY,Y, = 2/1/3SZF2, (8) 
y = -2(l-X){C^-X)^2Y^{C'-X)^X[Y^-CYl^S'Wl] 

= -2a-X){C^-X)-^2Y^{a-X)-^XY^[l-CHl-i-S'P,], (9) 

8 = -2aXSY,Ya^ = -2UsSaXY^, (10) 
and 

e = (1 -X){c^-x^f - ar\c^-x)-r-xs^,Yi 

= (l-X){C^-X)^-C^Y'iC^-X)-OSillXr. (11) 

Similar equations have been derived by [Johler and Walters, 1960] without neglecting the loss 
mechanism. 

For the present case I2 is set equal to zero to maintain the refracted angle in the plane of in- 
cidence. This is not an essential step but simply reduces the complexity of the solution. The 
more general case of /2 # would follow the procedures given without any complication. 

Once, a, /3, y, 8, and e are obtained, Booker's quartic equation, (6), may be factored and the 
four values of q obtained. The two roots, qi and ^2, which make the signs of ray path angles Og^ 
and Og^ the same as di are the ones being sought, since they represent waves travehng in the posi- 
tive 2-direction. 
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The angle Og is given by 

d§^ dy d8 d€ 

*"" '^ = Ml) ^ ^ 

These angles are the directions of two ray paths associated with q\ and ^2, respectively. 

The angles Op^ and dp^ of the refractive index n will also be useful in the following calculations. 
They are given by 

5 
tan^p--. (13) 

using the respective values q\ and 92- 

The relative phase of the ray at any point r is given by 

(t> = Sx-\-qz=n' r=nr cos (dg — Op), (14) 

and 

sin Op 

It is necessary to determine this relative phase since the phasor sum of the various scattered field 
components is to be computed to find the total scattered field. 

To specify the magnitude of the ray at the intersection point on the second interface, it is 
necessary to compute the transmission coefficients at the first interface. These are given by (A-1) 
and (A-2). One may also compute the reflection coefficient at the first interface at this time which 
would be associated with the directly reflected ray. 

The fields associated with the two rays at their point of intersection with the second interface 
(shown in fig. 2) are given by 

^UA) = Uo\\Tie-Mmr. cos iS^^-O^^)] (1^) 

and 

^r2A) ^ Uo\\T2e-J^^''''' ^^' (V-^^2)]' 

where n and r-z are the lengths of the ray paths. The directly reflected fields at the first interface 
are given by 

l|f/«(A)=f/„||R||oe-^w-t'-) (17) 



l|f/«(A) = f/o||«^oe--''*'^-^-''^'. 

It is now necessary to determine a new set of parameters associated with the transmission 
and reflection mechanism at the second interface. The first step is to transform the coordinate 
system such that the z[ axis is now normal to this second interface, as illustrated by x'., y' z'., 
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Figure 2. Ray tracing for nonparallel interfaces. 



coordinate system for ray 1 and shown in figure 2. The values q[ and S[ of the ray incident on 
this second interface in this primed coordinate system are [Lee, Peters, and Walters, 1964 a and b] 



cos [dp -h 0) 

f 

cos Or, 



(18) 



and 



S' = S 



sin {Op + 0) 
sin Op 



(19) 



The subscript (1) on the coordinate axes refers to the system associated with qi. A second coordi- 
nate system, designated by ^2, yg, ^2, is associated with qz for ray 2. This is necessary because 
the phase angle Op differs for the two ray paths n and r2. 

The value Ci^= 1 — S/^ and the values of the direction cosines // and W of the magnetic field 
are readily obtained in the new coordinate system. The coefficients of Booker's quartic equation 
(a', /3', y', 8', e') now may be obtained from (7-11). 

For ray 1, these coefficients are a function of S{ at point 0[. Thus Booker's equation must be 
factored to obtain the roots q[, q{2)[, q'^, q{A)[- However q[ is already known and Booker's equation 
reduces to the cubic equation 



a\q'-^ + (/S; + a\q'M'' + (y? + P[q[ + a[q['^)q' + (b[ + y[q[ + p\q\^ + a\q'-^) = 0, (20) 

where a[, j8J, y\, b[ are associated with S[. There are three solutions to (20) which are designated 
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by 9(2)1, ^3, and ^(4)1. It has been shown [Lee, 1963] that the desired value of ^3 should be chosen 
so that if q'i> q{2)'i, then q'^<qi4)'i, when the associated pair (^J, ^3) is being sought. Now (for 
ray 1) dg.^, shown in figure 2, is found from (12), and 6p'^ from (13). The index of refraction, n^, 
is found from (15). The same technique can be applied to find q'^ from q'^ at point O2 for ray 2. 

The point of intersection of the ray path with the next interface may then be found for the 
geometry being treated. In the example of figure 2, this is the original interface. Calculation of 
phase at this intersection follows the method described above. The reflection and transmission 
coefficients at this boundary, i.e., at the point O' of figure 2 are given by (A-3) and (A-4). At 
point 0;, q[, ^3, and q(4)[ are used to calculate [R'^, [T\\, and JT^. At point O'^, q'^, q'^, and ^(3)2 are 
used to calculate 2^4, '^\\^ and IT^. 

Thus only reflection coefficients [Ri and 2/^4 are considered in the following treatment. For 
cases where [R{4)[ and 2^(3)2 are not neghgible the additional computation [Lee, 1963] would parallel 
those used in this paper. Reflections described by coefficients [R{4)'i and 2^(3)2 are interpreted 
as coupling between waves. 

The fields transmitted into free space at the second interface are given by 

iiC^ii(i/.)==f/oiir, ;riie-^'^ie-^wi-i^Qn), ^21) 

l|f^ll('2z.) = f^o||r2 ^r||e--^'^^e-J-'^(^^'^'2-^c'^^^, 
and 

where 

(f)i = kTi\r\ cos{6p, — dg^); (\)2 = knir2 cosiOp^ — OgJ 

and all phases are referenced to the original origin of coordinates at point 0, i.e., x = y = z = 0. 
The fields of the ray reflected back to the gyrotropic medium to the point 0^ and O^ are 

UU = Uo\\T[[Rse-J^^r^^'^^ (22) 

and 

U^2B) = Uo\\T22R4e-J^^^^^'^\ 

where 

<i>[ = kn',r', cos {e;.-9-J and ct>'2 = kny, cos (01-0',), 

and all phases are referenced to the original origin 0. 

The use of ray optics for the gyrotropic body has considered a ray incident upon the body from 
the external medium and followed it completely through one internal reflection. AU of the fields 
associated with this case have been given. Any additional internal reflections may be treated 
simply by repeating these same steps. In addition, any coupling terms may also be readily in- 
cluded in any case where coupling becomes significant. 

Any changes in amplitude introduced by diverging ray systems, i.e., spatial attenuation, have 
been neglected. However, this problem can be handled by the modified geometrical optics method. 
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3. Conclusions 

A modified geometrical optics method has been described for determining the scattered fields 
of gyrotropic bodies, lliis method makes use of a ray optics technique developed for gyrotropic 
media. The general methods given are applicable to any anisotropic body. Other known tech- 
niques cannot be applied to find these scattered fields. 

4. Appendix. Reflection and Transmission Coefficients at Planar ''Free-Space 

Gyrotropic" Interfaces 

The coordinate system used in this appendix is chosen so that S\=S and S2 = 0. There are 
four possible w^aves in the gyrotropic medium; paired as upgoing waves, designated by subscripts 
(1) and (2), and downgoing waves, designated by subscripts (3) and (4). Thus a wave (1) or (2) 
in the gyrotropic medium may be reflected as wave (3) or (4), the reflection coefficient of which 
is designated i/?3, i/?4, 2R3, or 2R4, etc. 

4.1. Wave Incident From Free Space Up to Gyrotropic Medium 
For parallel polarized plane wave incidence {E'^ = 0) 

E' E' 9r 1 ^"^"^^ 

"^' i)H\^ ^ E[. Dr"""' C"^^'^' 
H'^ £•/; 1 1 



— C {7TyJ]y^ — r]y,TTy.,) " {TT j;TT y., " 7Ty,7Tj:J ] , 



and 



where 



£'f 2 

and the numerical subscript is the subscript for the particular q involved. 
For perpendicular polarized plane wave incidence (El — 0) , 

E^ 2 



T2--fr=^^i^^ACy]yX , (A-2) 
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-^^-' ='Fr = 7T" [~ ('^^.'^«2~'^J/.'^-<-2) +?; {■^X,VX^-VX,TTX2) 



— C{TryJ]y^ — 'r)y,7Ty^) + (TTx^TTy^ — TTy.TTx^) ] , 

and 






where 

Transmission into the gyrotropic medium yields coefficients of the form \\Ti^ ||7'2, where the 
subscript || designates the polarization of the incident wave and (1) or (2) designates one of the 
upgoing waves. Conventional definitions of polarization are used so that subscript || or x means 
that the E vector is "parallel to" or is "perpendicular to" the plane of incidence, respectively. 



4.2. Wave Incident From Gyrotropic Medium Up to Free Space 
For wave (1) (£'4 = 0) 

+ Vyi (-^x.Vx, - Vx,7Tx,H CTTx,{TTy,7)y, - 7]y,7Ty,) 

— CjTy, {'nx.J]y^ — Vvii'^xJ + C7]y, {7Tx,{TTy^ — 7Ty.^7TxJ } , 

E^ 1 



-^Vuii'^y:iVoc4—r}x:,'7Ty,) " "^ 7Tx^{7Ty^lqx^- Vx-.^TTyJ 



1 ( - \-l 



+ 77 TTyi {tTxJ]x^ — V-r-s'TTxJ ~-^Vxi {'^x-JTy^ — 7Ty,{TTx^)} •, 






^^'^^^r~^~Tr \-(vx,ny4-vyiVx,) +7^ (77^,17x4-^7^.7^^4) 



+ C {TTy,l)y, - y]y,TTy,) + {lTx,7Ty^ " 77^,77x4) ] , (A-3) 

and 

£/ 1 1 

1^4 = -^ = ^ [{VxrVy,-VyiVxs) +^ (77x117x3-^7x1^^3) 

-\-C{7Ty,7}y., — 7]y,7ryJ + (77j:, 77 j/3 " 77^,77 j,.,) ] . 
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Forwave(2)(£'/ = 0) 
r]H', 1 E''\ C 

+ V1I2 {tTx.Vx, — r)x:{TTx,) + CtTx^ {'7Ty,{Y]y, - r)y,,7TyJ 

— C7Ty^{TTx:{r]y^ — r]y.^7rxj + Cy]y.,{TTj:{Try, " TTy,{TTx,) }, 

E^ 1 

-^ '^yzi'^y^'^x^ — '^x.^TTy^) — — TTx^{TTy,{r]x., — 'r]x,7Ty^) 

+ Q '""^2 {^x.J]x^ — 'r)x:{TTx^ )—-^V-r2 ("TTx-.TTy,, — 7Ty,{TTx^ ) } , 
2R:i = ^ {'r]x.,'r)y,—'r]y.,nrx) -^Q {'^x.,y]x, — y]x.JTx,) + C{tT y.J]y, — y]y.JT y,) + {7rx.JTy,-TTy.JTx,) , (A-4) 



and 






^L z>2 



(V'-zVlh VuzVjJ +^ (^J^2T7^3 — T7./-277'.r:,) + C(77-,/,T7,/, — r)y.,7T y,,) + (Tr.r.TT,,, — 7r,/,7r.r,) , 



where 



Z>2 = {Vr-iVUi ~ Vy^Vx4) + 7^ (^^■a'»7x4 "" Vx.^^x,) -^C(7Ty,,7]y, — r]y,,7TyJ -\- {7Tx,7Ty, — TTy.JTx,) . 

The dominant reflection coefficients are i/?3 and iRa- Coefficients \Ra and 2/^3 are usually negligible. 



4.3. Wave Incident From Gyrotropic Medium Down to Free Space 



For wave (3) (^^=0) 



i]H'[ El -C 

3^11 = -^ = -C Y^=~fr {^-^,(^7 ./•.'^.'/2 - T7.V,T7x2) - 'r]x-,{'TTx,r]y., - ViuTTx.,) 



E' 



+ VU-k^J-xV-i-^ - VxyTTx^) — C7rx.,{Try,7]y., — r]y,7Ty.,) + CtT y.^TT x ,'r]y^ — 'r]y,'TTx^) " Cj]y,,[TTx,7T y., — TTy.TTx^)}, 



+ VyX'^yiVj^2 ~ Vx.TTy^) + ^ 7TxS7TyJ)x^ — y)x,7Ty^ 



~~F 7Ty^{7Txr'r)x2 ~ VxJTx^-\-j; '^x.i'TTx.'TTy^ — TTy.TTx.,) K 



E« 1 r 1 

^^^~E^^W \ ^^•^^'^^•■^ ~ Vy2Vx:) -Q ij^x2r)x^ - Vj^2^x.s) 



- C{TTyJ]y,, — 7)y.,7ry,,) + (TT x-JT y,, " TT y.JT x) L (A-5) 
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E^^ I { 1 1 



For wave {4^){EL=0) 



^a 



4 



inHl ^ El C ^ ^ ^ ^ X , / X 

— CTTxlTTyJ]y^ — riy.TTy^) + CTTySTTxyVy^ — Vy.TTx^ " Cr)y^{7Tx,7Ty, — TTyJTx^)}, 

EL. 1 

7 

2:4 



^^^^ff'^O'] '^ySy)xr'r]y2 - Vy^Vx^) — '^xSjTy{r]y^ — -qy.TTy^ 



+ Vy4i'^yiVj02 — Vxr^TyJ + ^ TrxJjTyJ]x^ — y]x,7Ty^ 

~~7^ '^yi'^xj]x2 ~ y]xJTx^^-?^ VjcS^ooi^Ty^ — TTy.TTx^ 

E^ \ { 1 

— CijTy^iqy, — iqy^TTyJ + {iTx^TTy^ — TTy^TTx^U (A-6) 

and 

4/?2 = -gT^ = ^ j - {rix,y]y, - VyiVj^J + ^ (^x,t7j-4 " t?^.^^^ + C{7ry,riy, — r)y,7Ty,) - (TTx^TTy, - ny.TTxJh 



where 



D3 = {rix,r)y2 ~ VyiV-r2) — Q ('^^•,^7^2 — '^X.TTx^ — CilTy.rjy^ - lf]y,7Ty^ + {TTx^TTy^ - TTy.TTx^)' 



In this case the dominant reflection coefficients are 3R1 and 4R2' The coefficients 3R2, 4R1 are 
usually negligible. 
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